Metastable States of a Coupled Pair on a Repulsive Barrier 
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cles with masses mi = mi = m, and coordinates ri and 
r-2 coupled by an oscillatory interaction through the po- 
tential barrier Vq{x{) + Vq{x2). The Hamiltonian of this 
system (% = 1) 
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Resonance penetration of two coupled particles through a repulsive barrier is considered. It 
is shown that a local minimum of the total potential generates metastable bound states, and their 
spectrum determines the position of resonances in the penetration probability. It is pointed out that 
the probabilities of tunneling of two interacting particles from the false vacuum can be essentially 
higher than it has been assumed earlier. 

PACS number(s): 03.65.Nk, ll.10.Jj, 21.45.+V 

In paper by N.Saito and Y.Kayanuma Q, it was 
pointed out that there exists a new quantum phenomenon 
- resonance transparency of a single repulsive barrier for 
a coupled pair of particles. To consider this effect, a 
one-dimensional rectangular repulsive barrier and an infi- 
nite one-dimensional rectangular potential well coupling 
the pair were chosen. Since the interactions were sim- 
ple, it was possible to solve the initial two-dimensional 
Schrocdinger equation by reducing it to a system of one- 
dimensional equations by means of projection onto 7 
eigenfunctions. However, there still remained the ques- 
tion of how this effect manifests itself in other systems. 

This note continues the study of the effect of reso- 
nance transparency of two type one-dimensional barri- 
ers for a pair of identical particles coupled by the os- 
cillatory interaction. This pair interaction allows us to 
reduce the problem of three-dimensional scattering of 
a three-dimensional oscillator to the solution of a two- 
dimensional equation analogous to the equation derived 
in ref. |lj]. Besides, it is just this sort of pair interac- 
tion that is used in literature || devoted to the proba- 
bility of induced decay of the false vacuum in collisions 
of high-energy particles (see, for instance, P,3|)- It was 
pointed out there that it is possible to describe the pro- 
cesses of transition from the false vacuum on the basis 
of quantum-mechanical tunneling of a pair of particles 
through the barrier; but the study was performed for a 
system where only one of the oscillator particles inter- 
acts with the barrier. Here, we show that, when the two 
particles interact with the barrier, there arises the same 
effect of resonance transparency as in ref. M . 

The first potential barrier we study is taken in the 
Gauss form from ref. || in order to show that it is possible 
to drastically increase the probability of induced decay 
of the false vacuum. The second potential barrier of the 
Coulomb form is investigated in order to draw attention 
to the fact that the resonance tunneling of the barrier is 
feasible in the problems of fusion of heavy nuclei. The 
method of investigation is based on the numerical solu- 
tion of the two-dimensional Schroedinger equation with- 
out any further simplifications. 

Consider the penetration of a pair of identical parti- 



written in coordinates of the center of inertia of the pair 
R = (ri + r2)/2 and in an internal coordinate of the rel- 
ative motion r = — r 2 describes the three-dimensional 
motion of a three-dimensional oscillator with the fre- 
quency of vibrations to. Since the potential barrier de- 
pends only on one variable, and the oscillatory interac- 
tion is additive in r projections, the wave function is fac- 
torized, and its nontrivial part describing the process of 
scattering depends only on two variables. It is convenient 
to represent these variables in the dimensionless form 



rmo . 



X2) 



— (X-L +X 2 ). 



The Schroedinger equation in these variables is of the 
form 

{-d 2 x -d 2 + x 2 + V(x -y)+V(x + y)-E)V = 0, 

(1) 

where the energy E is written in units u>/2, and the po- 
tential barrier V(x ± y) = ~Vo((x ± y)/\f2mto) in what 
follows will be written in a form convenient for us. Equa- 
tion (Q) should be supplemented with boundary condi- 
tions. Let the process of scattering proceed from left to 
right, and the initial state of the oscillator be the state 
n. Then the boundary conditions are written in the form 



exp(ik n y)(p n (x) - J2 S nj exp(-ikjy)ipj(x); 

Rnj exp(ik 3 y)(pj(x); 

0. 
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The wave functions of the oscillator (fj(x) obey the 
Schroedinger equation 



1 T 



(-flg + x 2 - ei) cpi = 
with the energy Ej = 2j + 1 (j = 0, 1, 2, . 



(3) 



momenta 



£j , and the number N of the last open chan- 



nel (E — Em+i < 0). In what follows, we consider an 
oscillator composed of bosons, whose spectrum is con- 
venient to number from 1. So, hereafter, Ej = Aj — 3 

U 

The probabilities of penetration W%j and reflection Djj 
are defined as the ratio of the density of a penetrated or 
a reflected flux to the density of the flux of incident par- 
ticles: 



W tj = 
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It is obvious that J2j<N (Wij + Dij) = 1 



To determine the probabilities of penetration (reflec- 
tion) in the above way, it is necessary to solve the two- 
dimensional differential equation (|l]); its numerical so- 
lutions will be given below. The numerical solution 
was based on a three-diagonal approximation of second 
derivatives with constant steps in x and y: h x = 0.025, 
h y = 0.005, respectively. Finite dimensions |y ma x| = 12 
and |a; max | = 7 of the range of numerical calculations, 
at a given degree of discretization, provided accuracy to 
within the third decimal point in all the presented calcu- 
lations. 

As it was indicated above, in papers devoted to the in- 
duced decay of false vacuum, use was made of the model 
of quantum-mechanical tunnelling of a pair coupled by 
the oscillatory interaction through a barrier || . The case 
when only one particle interacts with the barrier was con- 
sidered. In the framework of the resonance tunnelling, 
we can expect that the picture of tunnelling would essen- 
tially change when the interaction of both the incident 
particles with the barrier is switched on. The barrier 
used in ref. [Q], upon being made dimensionless, is of the 
form 



V(X) 



— expi-g'X'/uj); 
9 u 



X = x± y, 



(4) 



where u is the oscillator frequency, and g 2 <C 1 is the 
model parameter of false vacuum. In rcf. [|, the de- 
pendence of the tunneling probability was calculated in 
the range of g 2 from 0.09 till 0.01 and at fixed frequency 
lu = 1/2. In the present calculations, the same value of u 
but greater values of g 2 are accepted. The reason is that 
there arise extremely (in the framework of numerical cal- 
culations) narrow resonances in the energy dependence 
of the tunneling probability of a coupled pair through 
the barrier. Therefore, in Figjl], we present the results of 
numerical calculations of equation(^) with barrier (^) at 
g 2 = 0.5, 0.3, 0.2 denoted by letters A, B, C, respectively. 




FIG. 1. Total probabilities of penetration through the bar- 
rier. The solid curve is for potential (^J); the dotted one, for 
potential 2V(y). The explanation is given in the text. 



In this Fig., we draw the probabilities of penetration of 
a coupled pair from the ground state into all the possible 
states, i.e. W = J2j<N Wij at different values of g 2 . A 
prominent resonance dependence of the tunneling proba- 
bility shows that there do exist barrier resonances under 
discussion. Note that the first resonance with decreasing 
g 2 shifts towards higher energies, but the quantity g 2 E r 
diminishes. 

When g 2 = 0.2, the probability of resonance penetra- 
tion is many times (~ 10 s ) as large as the probability 
of penetration in the nonresonance region (background) . 
Therefore, in Fig.[j], we present the results of calculations 
on the logarithmic scale. They demonstrate not only the 
indicated exceeding but also the coincidence of the back- 
ground part of the curve with the probability of penetra- 
tion of a structureless particle, i.e. with the solution of 
equation (|l|) for the barrier potential 2V(y). To estimate 
the contribution of narrow resonances to the probability 
of penetration of the particle flux distributed over energy, 
in Fig.Q, we plot the penetrated flux 



1(E) 



1 



E — E< 



W(E')dE', 



in the case when the incident flux is distributed uniformly 
from Eo till E; the quantity Eq = 5. It is seen that 
the main contribution to the probability of penetration 
comes from resonances. At E = 23, the difference from 
the background penetration amounts to 4 orders. 

Now, we describe a simple scheme of arising bar- 
rier metastable states that make the barrier transpar- 
ent. It is not difficult to verify that the potential energy 
U(x, y) = V(x + y) + V(x — y) + x 2 possesses a local 
minimum at y — (the center of mass in the middle of 
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the barrier) and at some values of x = ±xo. A maximum 
is at x = 0. 

This, there exist two potential "wells" separated by the 
barrier. Bound states of that system split into even and 
odd states. The magnitude of splitting is determined by 
the probability of penetration through the internal bar- 
rier. When 2V(x = 0, y = 0) ^> 1, this shift can be very 
small, and the spectrum of even states is determined by 
the spectrum of an isolated "well". In the first approx- 
imation, the position of resonances can be described by 
the oscillator spectrum of bound states at y = and 

x — Xq 

E nx n v =E Q + 2w x (1/2 + n x ) + 2w y (l/2 + n v ), (5) 

where n x and n v are oscillator quantum numbers; Eq = 
2V(xq); and frequencies u> x and oj y are determined by 
the second derivatives at the point of local minimum: 

w x = ^/d 2 U(x,y)/2 , uj v = ^/d 2 U(x,y)/2. 

For potential (@), it is not difficult to obtain the 
oscillator-model parameters 

E =cu(l + 2ln(2/uj))/g 2 ; 

wl=Ahx(2/ U ); (6) 
^s=^ - !• 

Let us compare the positions of resonances drawn in 
Fig.nl at g 2 = 0.2 with the results of calculation by formu- 
lae (||) and (Q) from which it is clear that there exist small 
but clearly seen satellite resonances. The position of the 
first resonance E r = 13.65 is well described by the oscilla- 
tor energy in the ground state E 00 = 13.92. The second 
group of resonances is generated by a single excitation 
of oscillators either along y or along x: Eqi = 18.18, 
Eio = 18.63. They are associated with the resonances 
at energies 17.24 and 17.74. The third group of reso- 
nances is generated by a double excitation E 02 = 22.45, 
En — 22.90, E 2 o = 23.34, and respectively, resonance 
energies are 20.58, 20.88, 21.72. So, the simple oscillator 
model of a metastable barrier state gives a correct qual- 
itative picture of the origin of resonances. Comparison 
with Fig.0 shows that the largest values of the tunnel- 
ing probability correspond to metastable states with a 
minimal excitation along the coordinate of the center of 
inertia. 

Decomposition around the point of equilibrium does 
not exhaust all possibilities of the oscillator model. 
Agreement between the resonance energy and the energy 
of a metastable state can be improved by a simple varia- 
tional procedure. To this end, we consider the position of 
a minimum xq and frequencies u x and u> v to be unknown 
quantities that are determined by the minimum of the 
average of the total Hamiltonian 

H = (4> x 4> y \ -d 2 x - d 2 +x 2 + V(x -y) + V(x + y)\c}> y (f> x ) 
over normalized eigenfunctions of the oscillator in the 



ground state: 




Varying H over xq , uj x and u> y , we can derive a system of 
three nonlinear equations to be not presented here in view 
of their being cumbersome. The of the three equations 
are solved analytically: xq = xo(u) x ,u> y ), u> y = y/ui x — 1. 
In this way, the function of one variable H = H(w y ) is to 
be estimated numerically; its minimum determines the 
variational estimate E val for the first resonance. Note 
that the variational connection of oj x and uj y is the same 
as in the case of decomposition (^). In Table]!], we present 
the comparison of E vav with positions of the first reso- 
nance E r drawn in Fig|l]. Agreement can be considered 
good in the framework of the above-indicated accuracy. 



TABLE I. Comparison of positions of the first resonance 
with the variational estimate 



9 2 




E r 


g 2 E r uj/2 


0.5 


7.649 


7.62 


1.30 


0.3 


10.416 


10.38 


0.779 


0.2 


13.680 


13.65 


0.683 



When j 2 < 1, the variational expressions get simplified 
and allow the following decomposition: 

EZ~E 0() + O(g 2 ). 

It coincides, with an accuracy up to 0{g 2 ), with the 
energy derived by a simple decomposition around the 
minimum of U(x,y). So, the estimate of the resonance 
spectrum made by formulae (|^) and (|^) is asymptotic as 
g 2 — > 0. In particular, when g 2 —> 0, we can indicate the 
limiting position of the first resonance in units of g 2 , i.e. 
the quantity g 2 Eui/2 used in ref. ||: 

g 2 Euj/2 -> uj 2 (1 + 21n(2/w))/2. 

At oj = 1/2, this energy tends to 0.472; this is shown in 
the fourth column of Table |. In ref. Q where the study 
was made of the penetration of a pair through the bar- 
rier, a smooth curve was obtained for the probability of 
penetration in the energy interval from 1.2 to 2. From the 
presented calculations it follows that, if the interaction 
of both particles with the potential barrier is taken into 
account, this curve becomes essentially nonmonotone. 

The barriers considered above are of the form of the 
Gauss function. For completeness, below we present the 
calculations for a barrier of the Coulomb shape cutoff 
both at short and long distances: 

{Q / A m j n : \X\ < A m j n 
Q/\X\ : X miD < \X\ < X max ; X = x ± y. 
Q/X max • I A j > A max 

(7) 
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The cutoff at short distances was introduced for mod- 
elling the nuclear Coulomb barrier in the framework of 
constraints imposed by the one-dimensional scattering. 
With this cutoff, the notion of " barrier height" is mean- 
ingful for the one-dimensional model of scattering. For a 
greater analogy, the barrier width at \X\ = X m i n should 
be small in spatial units of the problem, i.e., as compared 
with the mean-square dimension of the oscillator. The 
cutoff at long distances was introduced to make it possi- 
ble to use an asymptotics of the type (g). The quantity 
I max should be larger than 1 for imitating the barrier of 
small transparency. Here we took X m i n = 0.1 X max = 5. 
The quantity Q determines the energy height of the bar- 
rier. In Fig.g, we show the results of calculations for Q = 
2, 4, 10 denoted by letters A, B, and C, respectively. 
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FIG. 2. Probabilities of penetration through barriers of the 
Coulomb type. Explanations are given in the text 

In this clear picture of the resonance tunneling 



of a coupled pair is also observed. We do not present the 
analysis of the oscillator model for the position of reso- 
nances, because the chosen potential is essentially of the 
model character. We only mention that satellite reso- 
nances manifest themselves clearly, and energies of prin- 
cipal resonances are equidistant. 

The considered mechanism of the transparency of bar- 
riers for a coupled pair manifests itself for all the poten- 
tial barriers chosen for the investigation. As the reso- 
nance transparency was first observed jj| for barriers of 
the rectangular shape, and the coupling in a pair was 
of nonoscillator type, it could be assumed that the res- 
onance transparency of barriers for composite particles 
could be observed for a wide class of interactions. There- 
fore, the effects of quantum transparency could occur in 
various fields of physics. In particular, when the interac- 
tion of two particles with a barrier is taken into account, 
the picture of induced decay of the false vacuum changes 
essentially. 
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the importance of the study of the resonance trans- 
parency of Coulomb barriers. 
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